library(maptpx)
## Loading required package: slam
library(qtlcharts)
library(data.table)
We thinned the gtex data so that the number of reads per tissue sample and per gene is few and of the order of the single cell level data. In the Jaitin et al dataset, there were on average, 2540 reads on an average and in the full experiment with 4590 cells and 20190 genes, the toal number of reads recorded was 11658921. On the other hand, for the GTEx dataset, the corresponding number was way bigger. For the 16407 cis genes, there were 3.5e+11 reads and for each sample on an average, there were 41073880 reads, which is around \(10^4\) larger than the single cell experiment reads in the Jaitin dataset. So, we use two thinning thresholds, \(p=0.0001\) and \(p=0.00001\).
## Set screen size to height=700 x width=1000
##
Read 60.9% of 16407 rows
Read 16407 rows and 8556 (of 8556) columns from 0.262 GB file in 00:00:05
Thinning parameter again \(0.001\).
##
Read 60.9% of 16407 rows
Read 16407 rows and 8556 (of 8556) columns from 0.262 GB file in 00:00:04